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Abstract 

We consider a primary model of ac-driven Brownian motors, i.e., a classical particle placed in 
a spatial-time periodic potential and coupled to a heat bath. The effects of fluctuations and dis- 
sipations are studied by a time-dependent Fokker-Planck equation. The approach allows us to 
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map the original stochastic problem onto a system of ordinary linear algebraic equations. The 
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solution of the system provides complete information about ratchet transport, avoiding such dis- 
advantages of direct stochastic calculations as long transients and large statistical fluctuations. 
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■ The Fokker-Planck approach to dynamical ratchets is instructive and opens the space for further 
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\£) ■ generalizations. 
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I. INTRODUCTION 



When one starts to pan for gold, he/she places a gold bearing sand (or gravel) into a pan 
and then agitates it in a circular motion, back and forth, left and right. Finally, the sand 
goes away with the water and gold grains remain in the bottom of the pan. However, if one 
simply inclines the pan then all the content will go away. 

This is the illustration of the basic mechanism which lies behind the ratchet idea: in 
a driven system without a preferable direction of motion, transport velocities depend on 
characteristics of movable objects - charges, spins, masses, sizes, etoi^ 1 ^. The intensive 
development of this idea during the last decade brought new tools for a smart control of 
transport in different systems, ranging from mechanical engines^ to quantum devices^. 

The ratchet setup demands three key ingredients which are (i) nonlinearity, (ii) asym- 
metry (spatial and/or temporal), and (hi) fluctuating input force of zero mean. The non- 
linearity, in a sense of a nonlinear system response to an external force, is necessary since, 
otherwise, the system will produce a zero-mean output from a zero-mean input. The asym- 
metry is needed to violate the left/right symmetry of the response because transport in a 
fully symmetrical system is unbiased. The zero-mean fluctuating force should break thermo- 
dynamical equilibrium, which forbids appearance of a directed transport due to the Second 
Law of Thermodynamics^. 

The basic model of dynamical ratchets is an one-dimensional classical particle in a pe- 
riodic potential exposed to an ac-held^^ 1 ^. While the symmetry analysis of microscopic 
equations of motion^ enables one to formulate necessary conditions for the directed transport 
appearance^, both the current sign and strength depend on dynamical mechanisms. At the 
deterministic limit, when noise is absent, the evolution of a damped particle is governed 
by attractors, regular (limit cycles) or chaotic ones 14 . Transport properties are encoded in 
the characteristics of attractors, so when there is only one attractor in phase space, the 
dc-current is equal to a mean velocity of the attractor—. 

The situation changes drastically when noise starts to contribute to the dynamics. During 
the evolution, particle jumps out of an attractor, evolves outside of the attractor vicinity, 
lands back into the attractor, etc. In other words, the particle explores the whole phase space 
and the dynamics of the particle can not be described in terms of the attractor properties 
only. The situation becomes even more complicated when several attractors coexist in phase 
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spaced. 

The standard approach, based on the direct Langevin simulations of the microscopic 
equations of motion 1 -, demands a very long time in order to overcome transient effects and 
to produce the sufficient self-averaging over a phase space. The goal of this paper is to show 
that all these problems can be tackled by using the time-dependent Fokker-Planck equation 
(FPE)iLi£. Our approach allows us to reduce the stochastic problem to a system of ordinary 
linear algebraic equations, which can be solved by using standard numerical routines. The 
solution of the system of equations provides a full information on transport properties of 
the system. 

The outline of this article is as follows. In Sec. II we introduce the model and set up the 
problem within the Fokker-Planck frame. Then, we formulate the symmetries which need 
to be violated in order to get a nonzero dc-current. Section III discusses the method of the 
solution of ac-driven Fokker-Planck equations. In Sec. IV, we illustrate our approach with 
the so-called tilting ratchets- Sec. V contains some conclusions and perspectives. 



II. THE MODEL 



The one-dimensional classical dynamics of a particle (e.g. a cold atom placed into an 
optical lattice, or a colloidal microsphere in a magnetic bubble lattice^), of mass m and 
friction coefficient 7, exposed to an ac-driven periodic potential and coupled to a heat bath, 
is described by the equation: 

rax + 7X = g(x, t) + g(x, t) = — dU(x, t)/dx, (1) 

where U(x,t) is a potential and the force g{x,t) is time and space periodic, 

g{x + L,t)=g(x,t + T)=g{x,t). (2) 

The noise is modeled by a (^-correlated Gaussian white noise, (£(£)) = 0, ((,(t)(,(t f )) = 
2 / y D5(t — £'), where the noise intensity is D = kT. 

The state of the system ([T]) can be represented as a point in the three-dimensional phase 
space, (x, v,t), for m 7^ (underdamped regime) and in the two-dimensional phase space, 
(x, t), when m = (overdamped limit )A= The stationary asymptotic current is equal to 

J = lim J(t) = lim (3) 
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The statistical description of the system at m = 1 is provided by the the Fokker-Planck 
equatior*i£ii£ 

{dt + iL v ~iL [iv - 9{x > t)] - iD & ]p(x > v > t] = °' (4) 

where v = x. The respective FPE for the overdamped limit, when inertia is neglible, m — 0, 
reads^ 1 ^ 

7 P(x, t) = -[^9(x, t) - D^}P(x, t) . (5) 

The overdamped limit is the appropriate description for microdynamics at low Reynolds 
number, when a particle moves in an extremely viscous medial. 

The Fokker-Planck equations, (HJ) , (E]) , are linear, dissipative, and preserve the norm, 
J Pdxdv for Eq. (Tj0) and J Pdx for Eq.([5]) 18 . In addition, the equations possess discrete 
time and space translation symmetries, so that the operations x — > x + L and t — > t + T 
leave the equations invariant. For a given boundary condition and a fixed norm, any initial 
distribution, P(...,0), will converge to a single time-periodic attractor solution, P(...,t) = 
P(..., t+T). What are the correct spatially boundary condition for the ratchet problem given 
by Eqs.([I]l3])? It has been shown that the current can be calculated with the spatial-periodic 
solution of the type P(x, ...) = P(x + L, ...)&. 

The dc-components of the directed current (jHJ) in terms of the spatially-periodic attractor 
solution P are given by^: 



J = ( / v ■ P(x, v, t) ■ dv)T,L, m = 1 , (6) 

J — oo 

J = 1 - 1 (g{x,t)-P{x,t)) T , L , m = 0, (7) 

where (...)t = Jq ■■■dt and = ...dx. Without loss of generality, we set L = 2tt. 

Let us assume that Eq.(j4j) (or Eq.(jSJ)) is invariant under some transformation of the 
variables x and t, which does not affect the boundary conditions. Then the unique solution 
P is also invariant under the transformation. The strategy is now to identify symmetry 
operations which invert the sign of the current J ([6]) (or (J7])) and, at the same time, leave 
the corresponding FPE invariant. If such a transformation exists, the dc-current J will 
strictly vanish 20 . Sign changes of the current can be obtained by either inverting the spatial 
coordinate x, or time t (plus simultaneous velocity inversion, v — > —v, for the underdamped 
case) . 
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Below we list all such transformations together with the requirements the force g{x,t) 
has to fulfill: 

S\ -x+x', t ^t + t'; Si{g) -> -g , (8) 

5 2 : x -> x + x', t -> -t + t'; S 2 (#)^£, (if 7 = 0). (9) 

Here x' and f depend on the particular shape of g(x,t). 
There is an additional symmetry for the overdamped limit", 

5 3 : x -> x + x', t -> -t + t'; S 3 (g)^-g, (ifm = 0). (10) 

which does not follow from the FPE (jSJ). This is not a symmetry of the equation (due to 
the last r. h. s. term in Eq.flSJ), so-called diffusive term 17 ). It has been shown, nevertheless, 
that the symmetry S3 corresponds to a certain property of the asymptotic attractor solution 
P(x,t), such that the current (jjTJ) exactly vanishes when S3 holds 20 . 

By a proper choice of g(x,t) all relevant symmetries can be broken, and then one can 
expect the appearance of a non-zero dc current J (JMZj)- 



III. METHOD OF SOLUTION 



It is impossible to solve the FPE's (jH-|3j) analytically. Therefore we deal with the equations 
by using the Fourier expansion (over the x and t) for the overdamped case (jSJ), 

N,K 

P{x,t) = -== Pnk-e< nM \ (11) 



v '2ttT 

n,k=—N,—K 

plus the matrix continued fraction technique^ 7 - for the underdamped regime (j3J), 

/ / \ N,K S 

P(x^t) = ^L J2 Y, P n*s-e iinX+kUJt) Mv), (12) 

^ n,k=-N-K s=0 

where ip s {v) is the Hermite function of s-order— . Both the expansions (jl 1111 2D are truncated, 
so we can control the convergence and the precision of the solution P by variation of the 
parameters N, K and S. 

By inserting the expansions (11111121) into the corresponding FPEs (EH), we obtain the 
following systems of linear algebraic equations, 

(iku + Dn 2 )P lk +in^2 9h p {n-i){k-q) = 0, (13) 

l,q 



for the overdamped limit, and 



(iku + sjD)P nks + iy/Dn(^P nk ( s ~i) + V s + lP nk ( s +i)) - J *5~] gi q P(n-i)(k-q)(s-i) = 0, 

V l,q 

(14) 

for the underdamped regime. Here g n k is the Fourier element of the force function, g(x, t) = 
k 9nke t( - nx+hujt ^ ■ We note that the equation for the dc-element, P 00 in Eq.(jSJ) (P oo in 
Eq. (HI)), should be replaced by the equation P 00 = 1 (Pooo = !)• This corresponds to the 
normalization condition, J^ 77 P(x,t)dx = 1 (J^ w J^° oo P(x,v,t)dxdv = 1). 
The definition for the current takes the form 



j = y^-Pooi, (is) 

for the underdamped case, and 

J = ~7= Yl 9nk ' P -n-k, (16) 

v n,k 

for the overdamped limit. 

In order to use the standard routine we have to transform the original variables, P n k 
(overdamped limit) and P n ks (underdamped regime), into the single- index variable V z , z = 
1, Z. There is the one-to-one transformation, 



{n,k} ->z= l + (n + N)(2K+l) + k + K (17) 
Z=(2N + 1)(2K + 1), (18) 

for the overdamped limit, and 

{n, k, s} -> z = 1 + s(2N + 1)(2K + 1) + (n + N)(2K + 1) + k + K (19) 

Z= {2N+1)(2K+1)(S+1), (20) 

for the underdamped regime, which transforms the original two- and three-dimensional 
matrices, P n k and P n k s , into the vector-column V z . 

The corresponding system of equations for V z can be solved by using the standard numer- 
ical routines. For our calculations, we have used the routine F07ANF from Linear Algebra 
Package (LAPACK) written in Fortran 77-22. 
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IV. EXAMPLE: TILTING RATCHETS 



There is a plenty of choices for a driving potential U(x, t) 4 . Below we consider the case 
of a particle placed onto a stationary periodic potential and driven by an alternating tilting 

g{x,t) = -V'{x)+E(t). (21) 

With the simple potential V(x) = V (l — cos(a;)), the two-frequency driving, 

E{t) = E Y sm{ut) + E 2 sm{2ut + 9), (22) 

ensures that for E\ , E% ^ S% is always violated. The symmetry S% is broken for 9 7^ 0, ±7r. 
In addition, S3 does not hold at 9 7^ ±7r/2 20 . 

We start from the overdamped limit, m = and 7 = 1. For a given set of parameters, 
Vo = 2, £?i = 4.6, £?2 = —6, and u = 0.75, there is only one limit cycle in the phase space 
of the system (inset in Fig.l). The limit cycle is nontransporting, i. e. particle returns to 
the initial position after one period of external driving T, for the whole range 9 G [— 7r, 7r]. 
Therefore, there is no dc-transport in the deterministic limit, D — 0. 

When noise enters the game, it changes the dynamics: particle can leave the attractor 
and this leads to a finite current appearance (Fig.l). Therefore, thermal fluctuations play 
a constructive role here by providing a way for the manifestation of the asymmetry hidden 
in the ac-driving. Since the current depends nonmonotonously on the noise intensity D, 
there is some kind of the stochastic resonance effect^ 4 -^. For a weak noise, the dynamics 
is still localized near the attractor (Fig. la, dashed line), so that the current is faint. 
In the opposite high-temperature limit, when the noise starts to dominate the dynamics, 
the particle dynamics is "smeared" over the phase space. The relevant space-temporal 
correlations are suppressed by the noise and the current tends to decrease. There is a 
resonance temperature (near D w 0.27, Fig. lb), where the dc-current can be resonantly 
enhanced. Here the directed flow of particles is analogous to the flow of information in the 
case of stochastic resonance^ 4 ^. 
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FIG. 1: (color online) (a) Dependence of the current J ()16p on 9, for the system ()13j) (N = K = 60), 
at different temperatures, D = kT. The parameter are 7 = 1, Vo = 2, E\ = 4.6, E2 = —6, and 
lo = 0.75. Insets: the attractor of the corresponding deterministic system ([I]) for 9 = — 7r/2(top 
left); the dependence of the running current J{t) = x(t)/t on t at the point 9 = —n/2, for D = 0.2 
(thick line) and D = 2 (dashed dotted line). Dashed lines correspond to the current values J (|16p 
(bottom right); (b) More detailed dependence: the current J (|16|) as a function of 9 and D. 

Note, that the noise variation can not violate the basic symmetries: current is always 
absent at the point 9 = and 9 = ±7r. We have corroborate our analysis by the direct 
Langevin simulation (inset on Fig. 1) and found a perfect agreement with the numerical 
solutions of the corresponding FPE. 
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FIG. 2: (color online) The attractor solution P(x,t) of the FPE (J5J, for L> = 0.02 (a) and D = 2 
(b). The white line in panel (b) marks the corresponding deterministic attractor. Here 9 = —ir/2, 
while the other parameters are the same as in Fig.l. 

The overdamped limit, however, is not suitable for the description of all realistic sit- 
uations. For example, for the modeling of Josephson ratchets it is necessary to take into 
account the inertia term, mas^. The underdamped regime without periodic driving has been 
studied for a long time. To the best of our knowledge, there is only one paper where FPE 
for an underdamped system with ac-driving has been solved numerically^. 

The results for m = 1, Vq = 1, Ex = E 2 = 1.2, u = 1, 7 = 0.1, and 6 = are shown in 
Fig. 3. There is a single chaotic attractor for the deterministic limit, D = 0. In Fig. 3a the 
Poincare section of the attractor (the position of the particle in the phase space, [x(t), v(t)], 
taken at stroboscopical times t n = nT 1 -), is depicted. Since the relevant symmetry, Si, 
is violated by the presence of the second harmonics E 2 , the attractor generates a finite 
dc-current, J « -0.0288 (Fig.4). 

Although the attractor solution of the corresponding FPE for D = 0.1 clearly resembles 
its deterministic counterpart (Fig.3), it produces a much stronger current, J ~ —0.083, than 
in the deterministic limit. Here, again, thermal fluctuations play a constructive role in the 
process of dc-current rectification. 
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FIG. 3: (color online) The stroboscopic representation of the attractor solution P(x, v, 0) (|12p of 
the Eq. (|14p (N = K = 35, S = 25), with the Poincare section of the corresponding deterministic 
attractor (fTJ) (white dots). The parameters are m = 1, Vq = 1, E\ = Ei = 1.2, u = 1, 7 = 0.1, and 
6 = 0. 

The convergence of the quantity (fl5l with increasing of the total number of basic states, 
Z pOj) . is relatively slow (Fig.4). We have to take into account 71 spatial and 71 temporal 
harmonics in order to estimate accurately the ratchet current. This fact means that the fast 
temporal modes (chaotic dynamics on a short time scale, t c T) and short spatial modes 
(x c L), contribute essentially to the overall ratchet current. 

V. CONCLUSIONS 

We have shown that the complex issue of ac-driven ratchets in a thermal environment can 
be resolved by using the Fokker-Planck equation. Our approach maps the original problem 
onto a set of ordinary algebraic linear equations, which then can be solved by a standard 
numerical routine. 

The approach opens a new perspectives for further developments. For example, it allows a 
generalization for the case of colored noise with exponential correlation function^. The cor- 
responding one-dimensional "non-thermal" process can be embedded into a two-dimensional 
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FIG. 4: (color online) The running current J{t) = x{t)/t vs t for deterministic system ([I]) and the 
dependence of the current J (|15p on the number of basic states Z (|20p . The other parameters are 
the same as in Fig. 3. 

"thermal" dynamics by introducing the additional dynamical variable, s(t), which should re- 
place the noise term, £(£), in Eq.(pQ)2£. The next prominent perspective is the quantum limit 
of underdamped ac-driving ratchets. The master equation for the corresponding Wigner 
function can be represented as a FPE (|14p with an additional coupling between elements^. 
Finally, the absolute negative mobility in ac-driven inertial systems^ 1 ^ is another possible 
target for our method. 
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